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In this paper, we evaluated the quasinormal modes of electromagnetic perturbation in a 
Schwarzschild black hole surrounded by the static spherically symmetric quintessence by using 
the third-order WKB approximation when the quintessential state parameter w q in the range of 
— 1/3 < w q < 0. Due to the presence of quintessence, Maxwell field damps more slowly. And 
when at —1 < w q < —1/3, it is similar to the black hole solution in the ds/Ads spacetime. The 

■ appropriate boundary conditions need to be modified. 
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I. INTRODUCTION 

(N ' 

The study of QNMs of black hole has a long history [H-[l2j since they were first pointed out by Vishveshwara flll i n 
calculations of the scattering of gravitational waves by a black hole. Many techniques such as WKB approximation 
[2fj| . the "potential fit" [2l|. and the method of continued fractions ( Leaver .1 985 |22j) have been developed to calculate 
^ the QNMs of black holes since Chandrasekahar and Deweiler[23| had succeeded in finding some of the Schwarzschild 
QN frequencies using the technique of integration of the time independent wave equation. In this paper, we use the 
, third-order WKB approximation. 

Recently. V. V . Kiselev |24j has considered Einstein's field equations for a black hole surrounded by the static spheri- 
\ cally symmetric quintessential matter and obtained a new solution that depends on the state parameter w q of the 
Q" 1 quintessence. And Songbai Chen et al[25] evaluated the quasinormal frequencies of massless scalar field perturbation 
around the black hole which is surrounded by quintessence. The result shows that due to the presence of quintessence, 
the scalar field damps more rapidly. 
ILJ ' In this paper, we discuss electromagnetic perturbation. 
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5-H ' II. ELECTROMAGNETIC PERTURBATIONS 

a ■ 

_ _ i 

The metric for the spacetime of Schwarzschild black hole surrounded by the static spherically-symmetric 
quintessence is given by[25j 

ds2 = (1 - ^ - ^Ti)* 2 - (1 - ^ ~ ^r 1 * 2 - + (1) 



where M is the black hole mass, w q is the quintessential state parameter, c is the normalization factor related to 
p q = ~ § r 3(i+L q ) i arid p q is the density of quitenssence. We consider the evolution of a Maxwell field in this spacetime. 
And the evolution is given by Maxwell's equations: 

F.p - 0, F^v = A Utfl - A^ v (2) 
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The vector potential A M can be expanded in four-dimensional vector spherical harmonics(see Ref.[26]): 
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Where I is the angular quantum number, m is the azimuthal number. The first column has parity (— and the 
second (—1)' 



Submit © to ©,and define *i = (1 - 2M - 



we can get 



— $(r) + {uj 2 - V)$(r) = 



(4) 



where $(r) = a Zm is for parity (— and $(r) = r 2 /Z(Z + 1)(— iujh lm — df lm /dr) is for parity (—1)' ( for further 

details see Ref.JH) 

The potential V is given by 



V(r) = (1 



2M 



c + 



3lL?q + l 
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III. WKB APPROXIMATION AND QUASINORMAL MODES 

The equation for perturbations of a black hole can be reduced to a second order differential equation in the form 



— <S> + (uj 2 -V)§ = 



(6) 



The coordinate x is a "tortoise coordinate" r* which ranges from — oo at the horizon to +oo at spatial infinity(for 
asymptotically flat spacetime).And the appropriate boundary conditions defining QNMs are purely ingoing waves at 
the horizon and purely outgoing waves at infinity. 
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-» +oo (8) 
Only a discrete set of complex frequencies satisfy these boundary conditions. 

B. F. Schutz and C. M. Will[l4[ presented a new semianalytic technique which is called WKB approximation for 
determining the quasinormal modes of black hole, later, S. Iyer and C. M. Will[l5[ developed the method to the third 
order and R. A. Konoplva[l6j extended it to the sixth order. The accuracy of the WKB formula is better with a larger 
multipole number I and a smaller overtone n. The formula for the complex quasinormal frequencies a; (for the third 
order WKB method) is 
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FIG. 1: Variation of the effective potential for the maxwell field in the Schwarzshild black hole surrounded by quintessence 
with r for I = 5, c = O.Of and w q = -1/300, -1/6, -0.33. 



and 



r,=r,(r p ) 



(11) 



For quintessence, w 9 keeps in the range of — 1 < w q < 0. Looking into the metric of we can see that, when w q is 
at the range of —1/3 < w q < 0, the spacetime is asymptotically flat. We can calculate the QNMs in the situation 
by WKB approximation. And when —1 < w q < —1/3, it is not asymptotically flat, which is similar to the ds/Ads 
spacetime. 

Take M — 1, c = 0.01 and M = 1, c = for our calculation. And c = means there is no quintessence. Using the 
third-order WKB approximation, we can get the solutions as the table 1 and table 2 show, where I is the angular 
harmonic index, n is the overtone number, u> is the complex quasinormal frequencies, w q is the quintessential state 
parameter. 



The variation of the effective potential with r which is respective to the quintessential state w q parameter for fixed 
I = 5 and c = 0.01 is shown in Fig. 1. 

TABLE I: The quasinormal frequencies of electromagnetic perturbations in the Schwarzshild black hole without 
quintessence(c=0) . 
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TAB IE ILThe quasinormal frequencies of electromagnetic perturbations in the black hole surrounded by quintessence 
for I = 2, I = 3, I = 4, I = 5 and c = 0.01. 



3w q + 1 


uj(n = 0) 


uj(ti = 1) 


Lu(n = 2) 


ui(n = 3) 


0.99 


0.45483 - 0.09458* 


0.43364 - 0.28950* 


0.40030 - 0.49335*' 


0.35869 - 0.70207*' 


0.8 
0.6 


0.45432 - 0.09442*' 
0.45364- 0.09421*' 


0.43318 - 0.28900*' 
0.43257 - 0.28833*' 


0.39992 - 0.49249*' 
0.39942 - 0.49134*' 


0.35842 - 0.70083*' 
0.35804- 0.69919*' 


0.4 


0.45280 - 0.09394*' 


0.43181 - 0.28749* 


0.39878 - 0.48990*' 


0.35755-0.69714*' 


0.2 


0.45175 - 0.09360* 


0.43086 - 0.28645*' 


0.39796 - 0.48811* 


0.35691 - 0.69459*' 


0.01 


0.45051 - 0.09321*' 


0.42972 - 0.28523*' 


0.39698 - 0.48602*' 


0.35611 - 0.69162*' 



3w q + 1 


w(n = 0) 


uj(n = 1) 


w(n = 2) 


u(n = 3) 


u(n = 4) 


0.99 


0.65343- 0.09515*' 


0.63825 - 0.28833* 


0.61202 - 0.48758* 


0.57849 - 0.69203*' 


0.53889- 0.89969* 


0.8 


0.65268- 0.09498*' 


0.63753- 0.28783i 


0.61137- 0.48673* 


0.57793 - 0.69082*' 


0.53842- 0.89810*' 


0.6 


0.65169 - 0.09477*' 


0.63660- 0.28717* 


0.61052 - 0.48560* 


0.57718- 0.68921*' 


0.53779 - 0.89600*' 


0.4 


0.65046 - 0.09449*' 


0.63543- 0.28634i 


0.60944 - 0.48418*' 


0.57622 - 0.68719*' 


0.53698 - 0.89337*' 


0.2 


0.64893 - 0.09415*' 


0.63396- 0.28530i 


0.60809 - 0.48242* 


0.57501 - 0.68467*' 


0.53593- 0.89010*' 


0.01 


0.64711 - 0.09376*' 


0.63222 - 0.28409i 


0.60648 - 0.48036*' 


0.57354- 0.68175*' 


0.53464 - 0.88630*' 



Zw q + l 


lu{h = 0) 


uj(n = 1) 


u(n = 2) 


uj{n = 3) 


uj(n = 4) 


0.99 


0.84873 - 0.09538*' 


0.83691 - 0.28788*' 


0.81544 - 0.48454*' 


0.78696 - 0.68574*' 


0.75317-0.89049*' 


0.8 


0.84774 - 0.09522*' 


0.83595 - 0.28738*' 


0.81454 - 0.48370*' 


0.78613 - 0.68454*' 


0.75242 - 0.88893*' 


0.6 


0.84645 - 0.09500*' 


0.83471 - 0.28672*' 


0.81335 - 0.48258*' 


0.78504- 0.68295* 


0.75143- 0.88686*' 


0.4 


0.84485 - 0.09472*' 


0.83314- 0.28589*' 


0.81187- 0.48118*' 


0.78366 - 0.68095*' 


0.75017-0.88426*' 


0.2 


0.84284 - 0.09438*' 


0.83119- 0.28485* 


0.81002 - 0.47943*' 


0.78192 - 0.67847*' 


0.74857-0.88102*' 


0.01 


0.84047- 0.09398*' 


0.82888- 0.28365*' 


0.80781 - 0.47738*' 


0.77985 - 0.67557*' 


0.74665 - 0.87726*' 



3w q + 1 


lu{h = 0) 


uj(n = 1) 


w(n = 2) 


Lu(n = 3) 


Lu(n = 4) 


0.99 


1.04260- 0.09550*' 


1.03293- 0.28765*' 


1.01484- 0.48279*' 


0.99013- 0.68164*' 


0.96033- 0.88391*' 


0.8 


1.04138- 0.09533*' 


1.03174- 0.28716*' 


1.01370 - 0.48196*' 


0.98905 - 0.68045*' 


0.95932- 0.88236*' 


0.6 


1.03980- 0.09512*' 


1.03018- 0.28650*' 


1.01220 - 0.48084*' 


0.98763 - 0.67888*' 


0.95799- 0.88031*' 


0.4 


1.03782 - 0.09484*' 


1.02824- 0.28567*' 


1.01032 - 0.47945*' 


0.98585 - 0.67689*' 


0.95632 - 0.87773*' 


0.2 


1.03534- 0.09450*' 


1.02581 - 0.28463*' 


1.00798 - 0.47771*' 


0.98360 - 0.67443*' 


0.95420 - 0.87453*' 


0.01 


1.03243- 0.09410 


1.02294- 0.28343*' 


1.00520 - 0.47567*' 


0.98094- 0.67155*' 


0.95167- 0.87079*' 



IV. DISCUSSION AND CONCLUSION 



The quasinormal modes of a black hole present complex frequencies, the real part of which represents the actual 
frequencies of the oscillation and the imaginary part represents the damping. 

From Fig.l we see that the peak value of potential barrier gets lower as the absolute w q increases. The data of table I is 
obtained by using WKB method in Schwarzschild black hole without quintessence, and table II is under the situation 
that the quintessence exits. Explicitly, we plot the relationship between the real and imaginary parts of quasinormal 
frequencies with the variation of ^(for fixed c = 0.01), compared with the situation without quintessence. From 
fig. 2 we can find that for fixed c(unequal to 0) and I the absolute values of the real and imaginary parts decrease as 
the absolute value of the quintessence state parameter w q increases. It means that when the absolute value of w q is 
bigger, the oscillations damp more slowly. And the absolute values of the real and imaginary parts of quasinormal 
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FIG. 2: The relationship between the real and imaginary parts of quasinormal frequencies of the electromagnetic perturbations 
in the background of the black hole surrounded by quintessence for fixed c = 0.01 and no quintessence for c = 0. x and *, 
respectively, refer to the value with quintessence and without quintessence. 



modes with quintessence are smaller compared with those with no quintessence for given I and n. That is to say, due 
to the presence of quintessence, the oscillations of the Maxwell fields damp more slowly. 

In this paper, we only calculate the QNMs when w q is at the range of— 1/3 < w q < 0. And for w q at the range 
of — 1 < w q < —1/3, the spacetime is not asymptotically flat, which is similar to the ds/Ads spacetime. And the 
boundary condition should be modified. An appropriate boundary condition is there are incoming waves at the inner 
horizon and outgoing waves at the outer horizon. Define the "tortoise coordinate" r* as ^ = (1 — — ^J^+i 
The boundary condition can be written as 

$~e-**"\r->r innei . (12) 

*~e <wr ',r- (13) 
We notice that when r = ri nner and r outeri r* is — oo and oo respectively. That is to say actually after the variation 
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of "tortoise coordinate", the boundary conditions for the —1 < w q < —1/3 case and the —1/3 < w q < case are 
formally the same. 
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